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Spin dynamics and fluctnations in the streaming regime 
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Spin dynamics of two-dimensional electrons in moderate in-plane electric fields is studied theoret¬ 
ically. The streaming regime is considered, where each electron accelerates until reaching the optical 
phonon energy, then it emits an optical phonon, and a new period of acceleration starts. Spin-orbit 
interaction and elastic scattering result in anisotropic relaxation of electron spin polarization. The 
overall spin dynamics is described by a superposition of spin modes in the system. The relaxation 
time of the most long-living mode depends quasi-periodically on the inverse electric field. The spin 
modes can be conveniently revealed by means of spin noise spectroscopy. It is demonstrated that 
the spectrum of spin fluctuations consists of peaks with the low-frequency peak much narrower than 
satellite ones, and the widths of the peaks are determined by the decay times of the modes. 

PACS numbers: 72.25.Hg, 72.25.Rb, 72.70.-|-m, 73.63.Hs, 78.47.db 


I. INTRODUCTION 

A moderately strong electric field applied to a semicon¬ 
ductor system with low carrier concentration can provide 
the streaming regime of electron transport [1]. In this 
regime, each free charge carrier accelerates quasiballisti- 
cally in the passive region of the momentum space, where 
its energy is smaller than the optical-phonon energy hujQ. 
As soon as the carrier energy amounts to it emits 
an optical phonon and scatters to a state with a small 
energy. This is the end of the period, after which the 
next cycle of acceleration starts. The electron momen¬ 
tum changes periodically from zero to a value po corre¬ 
sponding to the electron energy equal to fujo, Fig. 1. The 
period of oscillations is the electron travelling time from 

P = 0 to p = po: 

itr=Po/|ef|, (1) 

where e < 0 is the electron charge and £ is the electric 
field strength. The electron distribution in the momen¬ 
tum space is strongly anisotropic, i.e. its length in the 
field direction, po, is much larger than its width. The 
corresponding region of the momentum space is called a 
needle and is shown by the red solid line in Fig. 1. 

This regime of carrier transport has been studied in 
detail in three-dimensional systems. A semiconductor 
laser in p-Ge has been realized based on the streaming 
effect [2]. Many interesting features can be revealed by 
studying the electric current fluctuations in the stream¬ 
ing regime [3, 4]. Recent theoretical investigations of 
streaming in two-dimensional systems demonstrate col¬ 
lective wave-like excitations of the electrons with multi¬ 
branch spectra and considerable spatial dispersion [5]. 
Ballistic transport with dominant optical phonon scat¬ 
tering has been realized in graphene [6], and a number 
of interesting theoretical proposals have been made for 
graphene [7-9]. 

The inclusion of electron spin degree of freedom into 
the two-dimensional streaming-regime kinetics gives rise 
to rich spin-related phenomena. Due to electron drift in 
the electric field and linear in momentum Rashba and 



FIG. 1: Illustration of the streaming regime in the momentum 
space. During a ballistic motion inside the needle (red solid 
line), an electron spin rotates with the frequency fl{p) (blue 
arrows). The electron can be elastically scattered by an im¬ 
purity (star and dashed arrow), then it reaches the active area 
(red arrow), instantaneously emits an optical phonon (wavy 
arrow) and returns to p = 0 (dotted arrow). 


Dresselhaus spin-orbit interactions, electron spin process 
with an average frequency fldr- Electrically induced spin 
beats and long spin relaxation times as well as high de¬ 
gree of the current induced spin polarization have been 
predicted for such a system [10, 11] . However the compre¬ 
hensive study of spin dynamics in the regime of substan¬ 
tial spin rotations in each acceleration period (ildrt tr 1) 
has not been made yet. As we show, the spin dynamics is 
not reduced to a simple exponential relaxation, and the 
electron spin polarization at resonant conditions persists 
despite of multiple elastic scatterings. 
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Estimations show that the streaming regime in GaAs 
based heterostructures can be realized in reasonable fields 
f ~ 1 kV/cm, where the parameter fldr^tr is of the order 
of unity. This shows a possibility of experimental inves¬ 
tigations of spin-dependent phenomena in the streaming 
regime. 

The spin-dependent streaming kinetics has a few tem¬ 
poral ranges: ftr, the elastic scattering time Tp, and 
the spin relaxation time r^. Therefore it is natural to 
study spin dynamics in the frequency domain. This can 
be done by means of spin noise spectroscopy being a 
modern and very efficient tool for investigation of spin 
properties in various systems [12-14]. This method is 
based on the measurement of fluctuating spin signals in 
the ensemble of unpolarized carriers and allows simul¬ 
taneous resolving different time ranges [15]. The spin 
noise of free two-dimensional carriers in electric fields has 
been considered in nearly equilibrium conditions where 
the electric field is weak [16, 17]. Here we investigate spin 
fluctuations in heterostructures in the streaming regime 
which represents the opposite situation of strongly non¬ 
equilibrium electron gas. We demonstrate that the spec¬ 
trum of spin fluctuations in this case consists of a series 
of peaks with different widths. 

The paper is organized as follows. In Sec. II, the gen¬ 
eral theory of spin dynamics is developed for the two- 
dimensional streaming regime with account for spin-orbit 
interaction and elastic scattering. In Section III, we con¬ 
sider the spin dynamics in the presence of either longi¬ 
tudinal or transverse effective field. In Sec. IV, the spin 
noise in the streaming regime is investigated. Concluding 
remarks are given in Sec. V. 


II. SPIN DYNAMICS 


We address the streaming regime accounting for the 
spin-orbit interaction and elastic scattering by impuri¬ 
ties. We consider a semiconductor A 3 B 5 heterostructure 
grown along the axis z || [ 001 ] and choose the in-plane 
axes as x || [110] and y || [110]. In this coordinate frame, 
the linear in momentum Hamiltonian of spin-orbit inter¬ 
action takes the form [18] 

y-SO = P yx^yPx j^xy^xPy — 2 *^ • ( 2 ) 

where p is the two-dimensional electron momentum, ai 
{i = x,y) are the Pauli matrices, Pxy,Pyx are the spin 
splitting constants, and n(p) is the effective frequency 
of spin precession caused by Rashba and Dresselhaus ef¬ 
fects. 

Electron spin dynamics is described by a kinetic theory. 
The kinetic equation for the spin distribution function 


Sp{t) in the passive region reads [10, 11]: 




dt 


dpa 


= J dPp'Wpp'SpPt) + n{p) X Sp{t). ( 3 ) 


Here Wpp' = — p'^)/( 7 rrp) is the probability of elastic 

scattering from p' to p with Tp being the elastic scatter¬ 
ing time, and we assume the electric field to be applied 
opposite to the x || [lIO] axis. After each scattering, an 
electron accelerates in the x > 0 direction until reach¬ 
ing the border of the passive region, and then returns to 
p = 0. Such a trajectory is shown by arrows in Fig. 1. 

Spin relaxation in the streaming regime can be caused 
by both electron penetration into the active region 
{p > po) and by elastic scattering. The latter mechanism 
is shown to be somewhat more efficient [ 11 ], therefore, 
for the sake of simplicity, we will neglect the first mech¬ 
anism in this study. To that end we assume that the 
optical phonon emission time is infinitely short, so the 
spin distribution is nonzero only in the passive region 
{p < Po). Accordingly, electrons have zero energy and 
zero momentum immediately after optical phonon emis¬ 
sion. Hence we can separate the two contributions to the 
spin distribution function 

Spit) = S{py)9ipx)Slit) + s°^\p-1). (4) 

The first contribution describes the spin density in the 
needle {py = 0 ), while the second one stands for the spin 
distribution in all the passive region out of the needle. 

The streaming regime can be realized only if the elastic 
scattering is weak (r^ 3> ttr)i therefore we will consider 
spin dynamics up to the first order in the small parameter 
ttr/Tp <C 1. Moreover, we are aimed to solve the problem 
at the timescale ~ Tp or longer, therefore we assume that 
the majority of the carriers are in the needle at the time 
t — 0. Anyway this situation always establishes during 
the time ~ 2 ttr after the electric field is switched on. 
Hence the kinetic equations for the two components of 
the spin distribution read 




nip) X 


^ s;it) 

2'KpTp ■ 

(5b) 


Since the electron trajectories are closed in the p-space, 
the general solution of Eq. (5a) can be presented as a 
superposition of discrete spin modes: 


1 Ji) 

yn 


p 

n= — oo 1——1 ^ 


X exp 


- Pxttr /Po) - PxUr/ iPoTp) ■ ( 6 ) 
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Here n enumerates the modes, I = —1, 0,1 distinguishes 
different orientations of the normalized eigenvectors en\ 
uin are complex eigenfrequencies of the system, gn'^ are 
the coefficients, and 'R-{px) is the operator of rotation 
around the y axis by the angle 

Px 

‘^{pD = j '^^viPx) = Pyxttrpl/ifipo)- (7) 

0 

The appearance of the operator Tt{px) is caused by 
the fact that inside the needle the precession frequency 
^yiPx) is nonzero, and the electron spins are rotated in 
the (zx) plane. Once the spin distribution in the needle 

is known at t = 0, the coefficients gn'^ can be calculated 
as 

Po 

= J dpxe-'-'’i'^^^P-/P°ell^*TZ-\px)Sl{0), (8) 

0 


This condition reflects the fact of immediate optical 
phonon emission at reaching the border of the active re¬ 
gion (the right semicircle in Fig. 1). It is violated if the 
electron-phonon scattering is spin-dependent, but this ef¬ 
fect is estimated to be negligible in the streaming regime 
for typical structure parameters [10]. From the coupled 
set of Eqs. (9) and (11) we obtain: 

So(t)-S^^{t) = J ‘^5papSpit-ttr{pox-Px)/po), 

^ ( 12 ) 

where po = ipox,Py)- 

We solve this equation using the perturbation theory 
in the small parameter ftr/ip- la the absence of elas¬ 
tic scattering (ftr/^p = 0), it follows from Eq. (6) that 

the eigenfrequencies and eigenvectors en^ satisfy the 
equation 

7t(po)4'^=e-‘""’‘-e«. (13) 


where we have omitted the terms proportional to the first 
and higher powers of ftr/Tp- The coefficients gn'^ depend 
on excitation conditions: At resonant spin excitation in 
the vicinity of p = 0, gn'^ are of the same order for all n, 
while at non-resonant excitation the spin distribution at 
t = 0 is a smooth function of p, and gn^ drop with n as 
cx 1/n^. 

The spin distribution outside the needle, S°'^*{p;t), 
can be readily found by integration of Eq. (5b): 


5°"‘(p;t)= / dp'Jpp, 


Sp.{t-tt,{jPx-p'x)/po) 


POcc 




(9) 


Here 


POx — 


pI-p 


2 


and the tensor operator Qpp> denotes Green function of 
the ordinary differential equation 

where the Greek subscripts and superscripts denote the 
Cartesian components, and is the Levi-Civita sym¬ 
bol. The electrons are immediately taken off the left 
semicircle in Fig. 1 by the electric field, thus 

S°'^*{-pox,Py]t) = 0 . 


In order to satisfy this boundary condition we take the 
Green function, Gpp', to be zero at Px < p'x- 

In order to find the eigenfrequencies and spin modes 
in the streaming regime, one has to substitute general 
Eqs. (6) and (9) into the boundary condition 


Here 7^(po) is an operator of the electron spin rotation 
after the travel through all the needle. The corresponding 
rotation angle, see Eq. (7), is $(po) = ^^dr^tr with 


f^dr — (dyxP^i^ 

being the average precession frequency in the needle, 
hereafter we assume that jdyx >0. As a result, we find 
that the eigenfrequencies are combinations of multiples 
of the travel and drift frequencies: 

= I ( 27 rn/ttr - f^dr) , 

and Bn^ = ij, = (i ± \x)l\j2 with x,y,z being 

unit vectors along the Cartesian axes. Physically these 
frequencies reflect the fact that Sy{px) does not process, 
but oscillates in time with the period ftr- By contrast, the 
spin polarization in {zx) plane process with the frequency 
fldr in addition to the periodic oscillations. 

In the first order of the perturbation theory the eigen¬ 
frequencies become complex and can be presented as 
= ujn^ + Sn^ [19]. Provided 


f^drTp » 1 (15) 

to account for the elastic scattering, we substitute the 
eigenvectors Bn'^ into Eq. (6) find and then Eq. (12) 

yields the corrections SnK In the opposite case HdrTp < 1, 
the perturbation theory for degenerate levels should be 
used. 

The spin dynamics can be probed by Faraday, Kerr 
or ellipticity signals which are determined by the total 
electron spin, S(t). In the streaming regime most of the 
particles are localized inside the needle, therefore we have 


7 7 oo 1 

SSit)-S;^it)= J dpyS<’^\pox,Py;t). ( 11 ) sit)= dpxSlit)= Y. 

-PO 0 n=-oo; = -l 




(16) 
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Here we have introduced the average spin polarization in 
the nth mode 

Po 

(17) 

J Po 

0 

The results of this Section describe spin dynamics at 
arbitrary strong and anisotropic spin-orbit splitting. 

III. RESULTS AND DISCUSSION 

In general case the Green function Qpp> can not be 
found analytically, therefore in the next subsections we 
consider separately two limits: (i) the effective field (p) 
is oriented along a:-axis {Pyx = 0 , Pxy 7 ^ 0 ) and (ii) ll(p) || 
y {Pxy = ^,Pyx 7 ^ 0). In the end of this Section, we 
briefly analyze spin dynamics in the presence of both 
components in 0 (p). 

A. Longitudinal effective field f2(p) || x 

First we consider the limit of Pyx = 0, when the ef¬ 
fective field 11 (p) is parallel to the electric field. Since 
Hx cx Py, the spin precession in the needle is absent, and 
Eq. (13) is simplified to exp(—iwi^^ftr) = 1, which yields 

= 27rn/ftr- (18) 




Clearly in this limit the condition Eq. (15) fails, and be¬ 
low we apply the perturbation theory for degenerate lev¬ 
els. 

Out of the needle, electrons move ballistically con¬ 
serving the Py momentum component. Therefore, the 
electron spin rotates in the (yz) plane with the con¬ 
stant frequency 2PxyPy/h. Hence the Green function of 
Eq. (10), Gpp', is the identity operator for Sx, and it 
multiplies Sy{p;t) ±iSz{p-,t) in Eq. (12) by the factors 
exp [±iH 2 ;(pa; — p(,)/(e£)]. On average, scattered elec¬ 
trons have zero momentum py, therefore the spin polar¬ 
ization does not process even with account for scattering 
off the needle. Hence we can choose the basis vectors as 
=y,x,z for I = —1,0,1, respectively. Substituting 
these elP and cun^ from Eq. (18) into Eqs. (6) and (12), 
we obtain the corrections to the eigenfrequencies in the 
form 

i i /" d^p 

= -+ X- / — exp [ 27 rm (pox - Px + p) /po] 

Tp 2nTp J ppo 

X cos [2iPxyPy{p0x - Px)ttr/ifiPo)] ■ (19) 

Hereafter we assume that all integrations are performed 
over the passive region p < po only. 

Despite the developed formalism is quite cumbersome, 
its interpretation is straightforward. Let us consider the 
dynamics of x spin component. It corresponds to Z = 0, 


FIG. 2: (a) The spin distribution in the logarithmic scale 
for Tp = 3ttr, Ts = lOrp at the times t = 0 (black solid 
curve), 1.8rp (red dashed curve), lOr^ (blue dotted curve) 
and 20rp (green dash-dotted curve), (b) The decay rates of 
the first three 2 -spin modes for n = 0 (black solid curve), 
n = 1 (red dashed curve), and n = 2 (blue dotted curve) as 
functions of At/i, Eq. (21). The green dash-dotted line denote 
the spin splitting corresponding to the panel (a). The inset 
shows the decay rates of the x-modes (dots) and the analytical 
approximation Eq. (20) (solid line). 


and for a homogeneous distribution (n = 0 ) one finds 
= 0, i.e. the eigenfrequency is zero. This reflects 
the fact that the total spin polarization along the x axis 
is conserved, because it is not affected by the spin-orbit 
interaction. Therefore its dynamics is the same as for the 
particle distribution function. The decay rates of the ex¬ 
cited x-modes [Eq. (19) for / = 0 and n ^ 0] are nonzero. 
They are presented in the inset to Fig. 2(b) which shows 
that all of them are of the same order (~ 1/Tp). This 
means that the spin distribution relaxes during the time 
Tp to a constant (zero mode). In the limit of jnj 1 one 
can find from Eq. (19): 


- Im 



( 20 ) 





















5 


This analytical expression describes the decay of all 
modes except for n = 0 with accuracy of 4 %, see the 
inset to Fig. 2(b). 

The spin components S'®"* and S^y precess, which re¬ 
sults in spin relaxation. First we consider the limit where 
characteristic spin rotation angles out of the needle 


A(t) = PxyPottr/h ( 21 ) 


are small: Acj) <C 1. In this case all the excited modes 
relax with approximately the same rates as the a;-modes, 
i.e. ^ 1/Tp for n yf 0. The decay rates of the 

homogeneous spin distribution = —IhiSq^^^ are 

given by [20] 


1 _ 1 _ 7 (A^f 
r| 30 Tp 


( 22 ) 


B. Transverse effective field fl(p) || y 


If the effective field is perpendicular to the electric field, 
the dynamics of y spin component is the same as for the 
spin-independent distribution function. By contrast, the 
spin components lying in the [zx) plane feel the effective 
field fly oc Pa;, so the characteristic precession angles are 
not small even inside the needle where Px A Po- Therefore 
we concentrate on these components. 

At Pxy = 0, one can use the unperturbed basis vectors 
and eigenfrequencies and Cbn\- We find the Green 
function of Eq. (10) in the form Qppi = Tt{px)TZ~^{p'x), 
and from Eq. (12) we obtain 


^(0 



(23) 

Px + p) /po] ■ 


One can see that the spin relaxation time in this limit 
is much longer than Tp. This means that the spin distri¬ 
bution relaxes in two stages. First, after a time ^ Tp all 
the spin modes except for n = 0 decay, so the spin dis¬ 
tribution in the needle becomes nearly uniform. In the 
second stage, this uniform spin distribution decays with 
the rate 1 /t^. This two-stage relaxation is illustrated in 
Fig. 2(a). Interestingly, since the decay rates of all the 
excited modes are of the same order, the spin distribution 
almost conserves its shape during the relaxation, so that 
after the time ~ 2rp the initial shape of the distribution 
is still visible. 

The advantage of the presented theory is that it can 
describe spin dynamics and relaxation for arbitrary char¬ 
acteristic rotation angle Acj). If this parameter is of the 
order of unity, then all the spin modes decay with the 
rate ~ l/'T'p- The decay rates of the first three modes 
calculated by Eq. (19) are presented in Fig. 2(b). In the 
limit of Acj) —>■ oo, (5„ = —^jTp for any n, because a single 
elastic scattering is sufficient for complete spin dephas¬ 
ing. 

One can see that for the simple eigenfrequencies 
Eq. (18), the average spin polarization is nonzero only 
in the zeroth mode, since only the term with n = 0 
contributes to Sz{t), see Eq. (17). Therefore the total 
spin decays monoexponentially with the rate — Im(5g^^ 
for arbitrary values of j3xy despite in this case all the spin 
modes decay at different timescales. This spin relaxation 
rate is shown by the black solid curve in Fig. 2(b) as a 
function of the spin-orbit coupling strength. Note that 
the decay rates of all the modes are always smaller than 
Tp ^ because at least one scattering is needed to random¬ 
ize spin orientation. The series of local maxima in the 
spin relaxation rate presented in Fig. 2(b) roughly cor¬ 
responds to the condition Acj) = mr. Qualitatively, these 
maxima and the oscillations in the decay rates of the 
other modes, are related to the destructive interference 
of the spin rotation angles for scattered electrons. 


One can see that real parts of have opposite signs 
for I = ±1, while their imaginary parts coincide. 
Importantly, it follows from Eq. (23) that for all I 
6n\fldr) = ^o^^(Odr — 27rn/ttr) and, in particular, 

S^J^(2TTn/ttr) = 0 . 

This means that the nth spin mode does not decay at 
Odr = 27rn/ttr. On the other hand, it follows from Eq. (6) 
that, at Odr yf 0, the average spin polarization is nonzero 
in any mode. Therefore the decay time for the spin lying 
in the (zx) plane can be infinitely long. 

The fact that the spin relaxation time is infinite for 
Udi-itr = 27rn can be explained by considering an elec¬ 
tron’s motion from px to in the ballistic area. Its spin 
is rotated by the angle 4)(p(, ) — <l>(p^) around the y-axis, 
where $ ^ is given by Eq. (7). Since Py momentum 
component is constant at this motion, the rotation an¬ 
gle can be rewritten as 4)(p'^) — <l>(p^). Moreover, during 
spin-independent elastic scattering, p^ is conserved, and 
the spin direction is not changed. Therefore, if s is a spin 
of a given electron, then the quantity 

I={s,+ is;,)e-‘Wrtt.pVpg ( 24 ) 

is conserved during an electron’s motion in the passive 
region, and this is correct in all orders in ftr/Tp- Here we 
assume, as before, that the optical phonon emission time 
is infinitely short. Equation (24) demonstrates that when 
an electron reaches the active region, its spin is rotated 
around the y axis by the angle Hdr^tr irrespectively to a 
number of elastic scatterings. Therefore at Hdr^tr = 27rn 
the electron spin always returns to its initial direction 
after an optical phonon emission, i.e. spin relaxation is 
absent. This situation is similar to the persistent spin 
helix [21], but here it is realized in energy domain. The 
analogy comes from the fixed relation between the elec¬ 
tron displacement along the field direction, Aa:, and the 
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FIG. 3: (a) The decay rates of the first four spin modes, 

— Im((5l^Vp) for n = 0 (red solid curve), n = 1 (blue dashed 
curve), n = 2 (green dotted curve) and n = 3 (magenta dash- 
dotted curve). (b) Spin relaxation rate of the most long-living 
mode as a function of electric field (black solid curve) and its 
approximation Eq. (27) (red dashed curve). The spin relax¬ 
ation rate in the presence of both components of the effective 
field is shown by the blue dotted curve at = PyxI'i- 


gain of its energy /S.E = |e£|Aa; independent of the elec¬ 
tron trajectory. 

In general case the spin relaxation time, , corre¬ 
sponds to the smallest decay rate of all the modes: 

l/T-r = niin(- Im (5^1)) = - Im , (25) 

where 

n* = [r!drttr/(2^) + 1/2] (26) 

is the number of the most long living mode with [a;] being 
the integer part of x. The decay rates of the spin modes 
calculated after Eq. (23) are presented in Fig. 3(a). From 
this figure one can see that the spin relaxation time de¬ 
fined by Eq. (25) is a periodic function of Despite 

Ddr is the structure parameter, the travel time ftr can be 



FIG. 4: Average spin in the most long-living mode as 
a function of Ddrttr (black solid curve) and its asymptotic, 
Eq. (28) (blue dotted curve). Inset: Average spin in the nth 
mode, Sn, for n = 0 (red solid curve), 1 (blue dashed curve), 
2 (green dotted curve) and 3 (magenta dash-dotted curve). 


changed by the electric field. Figure 3(b) illustrates the 
oscillations of the spin relaxation rate as a function of 
the applied electric field in the streaming regime. In the 
limit Ddrttr ^ 1 we have from Eq. (23): 

4±i) « (2 + 4G - ^)-ii%^ (2^ + 1 - 6G), 

ZTTTp OTTTp 

(27) 

where G « 0.9159... is the Catalan constant. This 
asymptotic for the spin relaxation time is plotted in 
Fig. 3(b) by a red dashed curve. Comparison with the 
numerical calculation shows that this expression is cor¬ 
rect up to Ddi-ftr ~ 1- 

After decay of all the modes except for the n*th one, 
the total spin polarization is defined by two vector coef¬ 
ficients, see Eq. (16), and it lies in the (zx) plane. 

Their absolute values coincide and decrease as a function 
of n*. The inset to Fig. 4 shows the quantity 




,(i) 

'n, 


2 


as a function of the parameter Ddrttr for the first four 
modes. For Durttr = Oj average spin is present only in 
the zeroth mode. With increase of Ddrttr, the non-zero 
spin polarization arises in all the other modes due to a 
dependence of the spin precession frequency on p^- Fig¬ 
ure 4 shows js^Vj as a function of Ddrttr- One can see 
that the average spin polarization in the most long living 
mode exhibits damped oscillations. The jumps in the de¬ 
pendence js^Vj are related to the switching between the 
modes n and n -I- 1 at the points DdrGr = 27r(n -|- 1/2), 
see Eq. (26). For large Ddrttr the conserved spin decays 
as 



1 I TT 

2 y Ddrttr 


( 28 ) 
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This dependence is plotted in Fig. 4 by a dotted line. 

If both components, VLx and fly are nonzero, a numeri¬ 
cal solution of Eq. (10) is needed in order to describe the 
spin dynamics. This situation is qualitatively similar to 
the case of electric field orientation at arbitrary angle to 
the main axes. The results of calculations for /3xy ^ 0 
are shown in Fig. 3(b) by a blue dotted line. This de¬ 
pendence demonstrates that the nth mode decays even at 
f^dr^tr = 27rn if VLx is nonzero. Spin relaxation is switched 
on due to an additional phase oc I5xy acquired by electrons 
which violates invariance of I, Eq. (24). However, pro¬ 
nounced drops in the decay rate at resonant conditions 
f^dritr = 27rn are still present. 


IV. SPIN NOISE 

The analysis of the spin dynamics performed in the 
previous Sections reveals multiple timescales in the 
streaming regime: 


itr ^ '^s- 


Therefore studies in the frequency domain are useful, and 
the spin noise spectroscopy serves as an excellent tool to 
that end. In the steady state the majority of the electrons 
are in the needle, so the spin fluctuations are character¬ 
ized by the correlation functions 

+ (29) 

where angular brackets denote averaging over the time t 
for the given delay r. In the steady state, the distribu¬ 
tion function in the needle is uniform, so the one-time 
correlator is given by [22] 

(30) 

Here N is the electron two-dimensional density, and we 
ignore an electric field induced spin polarization because 
its value does not exceed a few per cents in the streaming 
regime [10, 11]. 

Ultimately we are interested in the total spin corre¬ 
lation function, {Sz{t)Szit + t)), where we assume the 
probe beam to propagate along z direction. Since a cor¬ 
relator of any physical quantity obeys the same linear 
kinetic equation as the quantity itself [3, 4, 22], the spin 
correlation function for r > 0 can be presented in the 
form 

{6Sz{t)6Sz{t + t)) 

= JdpxJ dp; J dp''(^6si^MT;.y,iT)6s;,,^^it)), 

where (r) {a, ji = x,y, z) is the Green function of the 
kinetic equation (3). Due to linearity of Tp^pi and using 


Eq. (30), this expression can be recast as 

{SSzit)SSzit + r)) = ^ ^SoSzir). 

p^p'^ 

(31) 

Here S'o = N/2, and Sz{t) is given by Eq. (16), where 
the spin distribution Sp^^z.ij) is found using the initial 
condition 


= So/po- (32) 

Accordingly the coefficients gn'^ can be calculated after 
Eq. (8): gn'’ = where is given by Eq. (17). 

The spin noise spectrum is defined by 

OO 

(SS^)^= J (5,(t)5,(t + r))e'“"dr. (33) 


Since the correlator is an even function of r [3, 4, 22], we 
obtain from Eqs. (31) and (33): 


N 


iss% = - 


4 - t 

n——(x> 1——1 




Im 


(0 

UJ — UJn 


U) ■ 


■ Wn 


(34) 

This equation demonstrates that the spectrum consists 
of the series of Lorentzian peaks centered at the eigenfre- 
quencies. The areas of the peaks are proportional to the 
squared total spin in the corresponding mode, and the 
widths are determined by the decay rates of the modes. 

At r2(p) II x, the average spin polarization along z- 
axis is nonzero only in the mode characterized by n = 0 
and 1 = 1. As it follows from Eq. (19), the correspond¬ 
ing eigenfrequency is imaginary, and = 1, see 
Eq. (17). Hence the spin noise spectrum simply reads 


{5Sl)^ = 


N 


2 1 + (t|w)^ 


(35) 


where l/r^ = —Im is shown by black solid line in 
Eig. 2(b). At 1, Tg is given by Eq. (22). 

In the opposite case of transverse effective field 
r2(p) II 2 /, all the spin modes with I = ±1 contribute to 
the spin noise spectrum. Figure 5 demonstrates the spin 
noise spectra for various values of Hdrttr- One can see 
that the spectrum consists of a series of peaks of differ¬ 
ent widths. 

If Hdrttr ^ 27r, the dominant contribution to the spin 
noise spectrum is given by the terms with n = 0 and 
I = ±1. Accordingly, at w > 0 the spin noise spectrum 
reads 




N t! 

4 1 -I- [(W - Hdr)Tj']^ ’ 


(36) 


i.e. it has a Lorentzian shape centered at the frequency 
Hdr- The width is l/r^ = — Im(5g^\ where is given 







defined in analogy with Eq. (33) has the form: 



FIG. 5: Spin noise spectra in the streaming regime calculated 
after Eq. (34) for Tp = 3ttr and = 0. 


by Eq. (27). The other peaks are centered at the fre¬ 
quencies ^1^71 jt^-^ ± fidr and have the widths of the order 
of l/r^; their amplitudes decrease as The shift 

of the main peak is caused by the nonzero effective pre¬ 
cession frequency inside the needle acting as a constant 
magnetic held. According to the fluctuation-dissipation 
theorem, this shift is equivalent to the electric-current 
induced shift of the electron spin resonance spectra [23]. 

For il<irttr/27r ~ 1, the widths of all the peaks have an 
order of I/t^, but the amplitudes of the hrst few peaks 
are comparable to each other, see Fig. 5. However, when 
Hdrttr approaches 27rn, the spin relaxation time tends to 
inhnity, see Sec. IIIB. As a result, the peak centered at 
w = |27rn/ttr ~ f2dr| becomes very high and narrow. 

Let us analyze the situation when the spin-orbit inter¬ 
action is absent. In this limit the total spin is conserved, 
and there are no spin fluctuations. Nevertheless, the spin 
noise spectroscopy has an access to spin dynamics even 
in this case. Indeed, in transmission or reflection experi¬ 
ments, the measured spin Faraday or Kerr signals, 0(t), 
are related to the spin distribution function via 


0(t)= I d^pK{p)SpAt)- 


(37) 


Here K{p) oc (p^ -|- a^)“^, where a is determined by the 
detuning between the probe beam frequency and the en¬ 
ergy gap [11]. The resonant dependence K{p) reflects 
the fact that the electrons with larger energy give smaller 
contribution to the spin signals, and it gives rise to fluc¬ 
tuations of 0(t). Analysis shows that the spectrum 0j 
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(“*)" = T E I® 


' Im 


n— — oo 
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where 

Po 

0„ = /^K(p)e2™P/^’“. (39) 

J Po 
0 

The eigenfrequencies cjn'^ = 27rn/ttr + where 
should be calculated after Eq. (19). One can see that the 
noise spectrum 0^ has a structure of Lorentzian peaks 
centered at multiples of the travel frequency and having 
the widths of the order of 1/rp (except for n = 0). This 
spectrum is different from the electric current fluctuation 
spectrum [3, 4, 19] by the zero frequency peak: since the 
spin relaxation time is infinite in this limit, the width of 
the peak is zero. The above analysis demonstrates that 
the proposed method extends the spin noise spectroscopy 
technique to the case when the total spin is conserved and 
does not fluctuate. The same approach allows measur¬ 
ing, e.g. energy relaxation rate of free electrons if spin 
relaxation is slow enough. 


V. CONCLUSIONS 

We have developed a kinetic theory of spin dynamics 
in the streaming regime with account for elastic scat¬ 
tering and spin-orbit interaction. The spin eigenmodes 
are identified and their decay rates are calculated. We 
have shown that electron spin dynamics is strongly dif¬ 
ferent in the limits of small and large spin rotation angle 
during the time ftr- If it is small, then the spin distribu¬ 
tion becomes uniform inside the needle on the timescale 
of one elastic scattering event. Afterwards, the average 
spin polarization monoexponentially decreases with the 
decay time Tg S> Tp. In the opposite limit of large rota¬ 
tion angles, the spin relaxation time has an order of Tp. 
However, the spin relaxation time oscillates as a function 
of the electric held and inhnitely increases when Hdrttr 
approaches a multiple of 27r. We have demonstrated that 
this effect is robust against elastic scattering, and, in fact, 
it is the energy space analogue of the persistent spin helix. 
The pronounced oscillations exist even in the presence of 
the transverse component of the effective magnetic held. 

The spin noise spectrum in the streaming regime is cal¬ 
culated. The spectrum consists of a series of the peaks 
corresponding to the different spin modes, and the widths 
of the peaks correspond to the lifetimes of the modes. 
The present study demonstrates that the spin noise spec¬ 
troscopy applied to nonequilibrium electron systems re¬ 
veals the parameters of spin dynamics and, in particular, 
the spin-orbit splittings. In the range of low frequencies 
inherent to the traditional spin noise spectroscopy, evo¬ 
lution of the spectrum rehects the strong oscillations of 
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the spin relaxation rate. The advantage of the ultrafast 
spin noise spectroscopy [24] paves the way for observation 
of peaks in the spin fluctuation spectrum in the stream¬ 
ing regime. Moreover, if the spin-orbit coupling is small, 
the resonant measurement of the Faraday or Kerr angle 
fluctuations allow investigating the particle distribution 
function dynamics. 

As an outlook we note, that spin dynamics in the 
streaming regime is extremely interesting to investigate 
in topological insulators. One of the reasons is a high 
value of the current-induced spin polarization. It is estab¬ 
lished that the spin polarization in topological insulators 
is proportional to a ratio of the drift and the Fermi mo¬ 
menta which is small in weak fields [25-28]. By contrast. 


in the streaming the needle-like electron distribution re¬ 
sults in a large drift momentum, and the current-induced 
spin polarization is 100 %. 


Acknowledgments 

We thank M. M. Glazov for fruitful discussions. Partial 
support from RFBR and RFBR-DFG ICRC TRR160, 
Dynasty Foundation, RF President Grant No. SP- 
643.2015.5 and Programmes of RAS is gratefully ac¬ 
knowledged. 


[1] A. A. Andronov, in Spectroscopy of nonequilibrium elec¬ 
trons and phonons, edited by C. V. Shank and B. P. 
Zakharchenya (Elsevier Science Publishers B.V., 1992), 
p. 169. 

[2] E. Gornik and A. A. Andronov, Special Issue on Far- 
Infrared Semiconductor Lasers, Opt. Quant. Electron. 
23 (1991). 

[3] V. Bareikis, R. Katilius, J. Pozhela, S. V. Gantsevich, 
and V. L. Gurevich, in Spectroscopy of nonequilibrium 
electrons and phonons, edited by C. V. Shank and B. P. 
Zakharchenya (Elsevier Science Publishers B.V., 1992), 
p. 327. 

[4] Sh. Kogan, Electronic noise and fluctuations in solids 
(Cambridge University Press, 2008). 

[5] V. V. Korotyeyev, V. A. Kochelap, and L. Varani, Wave 
excitations of drifting two-dimensional electron gas under 
strong inelastic scattering, J. Appl. Phys. 112, 083721 
( 2012 ). 

[6] A. Barreiro, M. Lazzeri, J. Moser, F. Mauri, and A. 
Bachtold, Transport properties of graphene in the high- 
current limit, Phys. Rev. Lett. 103, 076601 (2009). 

[7] T. Fang, A. Konar, H. Xing, and D. Jena, High-held 
transport in two-dimensional graphene, Phys. Rev. B 84, 
125450 (2011). 

[8] S. Sekwao and J.-P. Leburton, Hot-electron transient 
and terahertz oscillations in graphene, Phys. Rev. B 83, 
075418 (2011). 

[9] S. Sekwao and J.-P. Leburton, Soft parametric resonance 
for hot carriers in graphene, Phys. Rev. B 87, 155424 
(2013). 

[10] L. E. Golub and E. L. Ivchenko, Spin-dependent phe¬ 
nomena in semiconductors in strong electric helds. New 
J. Phys. 15, 125003 (2013). 

[11] L. E. Golub and E. L. Ivchenko, in Advances in Semi¬ 
conductor Research: Physics of Nanosystems, Spintron- 
ics and Technological Applications, edited by D. Per- 
sano Adorno and S. Pokutnyi, (Nova Science Publishers, 
2014), p. 93. 

[12] V. S. Zapasskii, Spin-noise spectroscopy: from proof 
of principle to applications, Adv. Opt. Photon. 5, 131 
(2013). 

[13] J. Hiibner, F. Berski, R. Dahbashi, and M. Oestreich, 
The rise of spin noise spectroscopy in semiconductors: 
From accoustic to GHz frequencies, Phys. Status Solidi 


B 251 , 1824 (2014). 

[14] V. S. Zapasskii, A. Greilich, S. A. Grooker, Yan Li, G. G. 
Kozlov, D. R. Yakovlev, D. Reuter, A. D. Wieck, and 
M. Bayer, Optical spectroscopy of spin noise, Phys. Rev. 
Lett. 110, 176601 (2013). 

[15] D. S. Smirnov and M. M. Glazov, Exciton spin noise in 
quantum wells, Phys. Rev. B 90, 085303 (2014). 

[16] F. Li, Y. V. Pershin, V. A. Slipko, and N. A. Sinitsyn, 
Nonequilibrium spin noise spectroscopy, Phys. Rev. Lett. 
Ill, 067201 (2013). 

[17] V. A. Slipko, N. A. Sinitsyn, and Y. V. Pershin, Hybrid 
spin noise spectroscopy and the spin hall effect, Phys. 
Rev. B 88, 201102 (2013). 

[18] N. S. Averkiev and L. E. Golub, Spin relaxation 
anisotropy: microscopic mechanisms for 2D systems, 
Semicond. Sci. Tech. 23 , 114002 (2008). 

[19] I. B. Levinson and A. Yu. Matulis, Gurrent huctuations 
in strong electric helds, Sov. Phys. JETP 27 , 786 (1968). 

[20] Equation (22) corrects a misprint in a factor of two made 
in Ref. [11]. 

[21] B. A. Bernevig, J. Orenstein, and S.-C. Zhang, Exact 
SU(2) symmetry and persistent spin helix in a spin-orbit 
coupled system, Phys. Rev. Lett. 97 , 236601 (2006). 

[22] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics 
(Butterworth-Heinemann, 2012). 

[23] Z. Wilamowski, H. Malissa, F. Schafller, and W. Jantsch, 
g-factor tuning and manipulation of spins by an electric 
current, Phys. Rev. Lett. 98, 187203 (2007). 

[24] F. Berski, H. Kuhn, J. G. Lonnemann, J. Hiibner, 
and M. Oestreich, Ultrahigh bandwidth spin noise spec¬ 
troscopy: Detection of large g-factor huctuations in 
highly-n-doped GaAs, Phys. Rev. Lett. Ill, 186602 
(2013). 

[25] L. E. Golub and E. L. Ivchenko, Spin orientation by elec¬ 
tric current in (110) quantum wells, Phys. Rev. B 84, 
115303 (2011). 

[26] T. Misawa, T. Yokoyama, and S. Murakami, Electromag¬ 
netic spin polarization on the surface of topological insu¬ 
lator, Phys. Rev. B 84, 165407 (2011). 

[27] D. Pesin and A. H. MacDonald, Spintronics and pseu- 
dospintronics in graphene and topological insulators, 
Nat. Mater. 11 , 409 (2012). 

[28] C. H. Li, O. M. J. van’t Erve, J. T. Robinson, Y. Liu, 
L. Li, and B. T. Jonker, Electrical detection of charge- 



10 


current-induced spin polarization due to spin-momentum 
locking in Bi 2 Se 3 , Nat. Nanotech. 9, 218 (2014). 


